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We use three dimensional time-dependent Schrodinger equation (3D-TDSE) to calculate angular 
electron momentum distributions and photoelectron spectra of atoms driven by spatially inhomo- 
geneous fields. An example for such inhomogeneous fields is the locally enhanced field induced 
by resonant plasmons, appearing at surfaces of metallic nanoparticles, nanotips and gold bow-tie 
shape nanostructures. Our studies show that the inhomogeneity of the laser electric field plays an 
important role in the above threshold ionization process in the tunneling regime, causing significant 
modifications to the electron momentum distributions and photoelectron spectra, while its effects 
in the multiphoton regime appear to be negligible. Indeed, through tunneling ATI process, one can 
obtain higher energy electrons as well as high degree of asymmetry in the momentum space map. 
In this study we consider near infrared laser fields with intensities in the mid- 10 14 W/cm 2 range 
and we use linear approximation to describe their spatial dependence. We show that in this case 
it is possible to drive electrons with energies in the near-keV regime. Furthermore, we study how 
the carrier envelope phase influences the emission of ATI photoelectrons for few-cycle pulses. Our 
quantum mechanical calculations are fully supported by their classical counterparts. 

PACS numbers: 42.65.Ky,78.67.Bf, 32.80.Rm 



I. INTRODUCTION 

The process known as above-threshold ionization 
(ATI), in which an atom or molecule absorbs more pho- 
tons than the minimum number required to single ionize 
it, has been a subject of intensive studies during the last 
decades (see e.g. [1] and references therein). The first 
experimental realization was made at the end of the sev- 
enties Q and since then there has been a truly amazing 
progress in understanding of the non-perturbative nature 
of ATI. The recent advances in laser technology made 
possible to routinely generate laser pulses with few cy- 
cles of duration, which allows control of the atomic and 
molecular processes in their natural time scales, i.e. in 
the range of (sub)-femto to attoseconds. In addition, 
these short laser sources find an extensive range of ap- 
plications in basic science, such as controlling molecular 
motions and chemical reactions Furthermore, the 

few-cycles pulses provide the fundamental pillar in the 
generation of high-order harmonics and the creation of 
isolated extreme ultraviolet (XUV) pulses 0] . 

The appearance of COLTRIMS experiments (see 
e.g. @ and references therein) offered an unprecedented 
possibility of performing stringent tests on the differ- 
ent theoretical approaches. On one side, this is because 
the imaging of the vectorial momentum distributions of 
the reaction fragments are easily accessible, while on the 
other they are particularly sensitive to various details of 
the theory. COLTRIMS were primarily developed for 
the study of few-body dynamics induced by particle im- 
pact, i.e. electrons and ions, but the extension to scruti- 
nize and tackle laser-induced processes was natural (see 
e.g (8l4To|). Among the features which were theoretically 



analyzed was the complex emission pattern present in the 
two-dimensional momentum plane, parallel and perpen- 
dicular to the laser polarization axis, of the laser- ionized 
electron distributions near threshold (TJ. It was also in- 
vestigated how these patterns evolve as the laser-matter 
process change from the multiphoton to the tunneling 
regimes [Tof . 

The main difference between a few-cycle pulse and 
a multicycle one is the strong dependence of the laser 
electric field on the so-called Carrier Envelope Phase 
(CEP) [H, [l3|. The electric field in a few-cycle pulse 
can be characterized by its duration and by the CEP. 
The influence of CEP has been experimentally observed 
in high-harmonic generation (HHG) [H| , the emission di- 
rection of electrons from atoms [l5|, and in the yield of 
nonsequential double ionization [16j. In order to have a 
better control of the system on an attosecond temporal 
scale it is, therefore, important to find reliable and di- 
rect schemes to measure the absolute phase of few-cycle 
pulses. 

The investigation of ATI generated by few-cycle driv- 
ing laser pulses plays a key role in the CEP character- 
ization due to the sensitivity of the energy and angle- 
resolved photoelectron spectra to the value of the laser 
electric field absolute phase [13, G3 . Consequently, the 
behavior of the laser-ionized electrons renders the ATI 
phenomenon a very valuable tool for laser pulse charac- 
terization. To determine the CEP of a few-cycle laser 
pulse, it is essential to record the difference between the 
yield of electrons ionized for different emission angles. 
Through this technique one can analysis the so-called 
backward-forward asymmetry in order to obtain the ab- 
solute CEP [l9|| . Furthermore, it appears that the high 
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energy region of the photoelectron spectra is most sen- 
sitive to the absolute CEP and consequently electrons 
with large kinetic energy are needed in order to describe 
it 0,|23 . 

Recent experiments using a combination of plasmonic 
nanostructures and rare gases have demonstrated that 
the harmonic cutoff of the gases could be extended fur- 
ther than in conventional situations by using the field, 
locally enhanced due to the coupling of a laser pulse with 
a nanosy stems [HI, [22|. In such nanosy stems, due to the 
strong confinement of the plasmonics spots and the dis- 
tortion of the electric field by the surface plasmons, the 
locally enhanced field is not spatially homogeneous in 
the region, where the electron dynamics take place. One 
should note, however, that the outcome of the experi- 
ments of Ref. [21} , in which a combination of gold bow- 
tie nanostructure and argon gas for generation HHG was 
used, has been recently under intense scrutiny [23|, l24j . 
In addition, recently, instead of atoms or molecules in gas 
phase solid state nanostructures have been employed as 
a target to study the photoelectron emission by few in- 
tense laser pulses [25, 26]. This laser driven phenomenon, 
called Above Threshold Photoemission (ATP), has re- 
ceived special attention due to the novelty of the involved 
physics and potential applications. In ATP process, the 
emitted electrons have e nergy far beyond the usual cutoff 
for noble gases (see e.g. j26l431j). Furthermore, the pho- 
toelectrons emitted from these nanosources are sensitive 
to the CEP, and consequently this fact plays an impor- 
tant role in the an gle and energy resolved photoelectron 
spectra 0, EaBl . 

From theoretical point of view, the fundamental as- 
sumption behind strong field phenomena, that the laser 
field is spatially homogeneous in the region where the 
electron dynamics takes place [34], [35| , is not any more 
valid for the locally enhanced plasmonic field. Indeed, 
in such system the driven electric field, and consequently 
the Lorentz force the electron feels, will also depend on 
position. Up to now, there have been very few studies in- 
vestigating the strong field phenomena in such spatially 
inhomogeneous fields |36i-l42j. All of these studies have 
demonstrated that the spatial dependency of the field 
strongly modifies the laser-driven phenomena that ap- 
pear in such circumstances. 

For homogeneous driving field, up to know, different 
numerical and analytical approaches have been employed 
to calculate the ATI (see e.g. p], l43l - [47| and references 
therein). In this article, we extend the studies of our pre- 
vious paper [5o| by applying the numerical solution of the 
time-dependent Schrodinger Equation (TDSE) in three 
dimensions to calculate the angular electron momentum 
distributions and photoelectron spectra of ATI driven by 
spatially inhomogeneous fields, considering both tunnel- 
ing and multiphoton regimes. The spatial dependence of 
the field is considered to be linear. We mainly focus on 
studying ATI of hydrogen atoms, but our scheme, within 
the single active electron approximation, can be directly 
applied to any complex atom. We demonstrate how the 



inhomogeneity of the field modifies both the two electron 
momentum distributions and the photoelectron spectra; 
we examine also the influence of the CEP parameter. 
Furthermore, our quantum mechanical results are com- 
pared with classical calculations of the kinetic energy of 
the electron. 

This article is organized as follows. In the next sec- 
tion, we present our theoretical approach to model ATI 
produced by spatially nonhomogeneous fields, with main 
emphasis on the extraction of the electron angular mo- 
mentum distributions starting from the TDSE outcomes. 
Subsequently, in Sec. Ill, we apply our method to com- 
pute the electron momentum distributions and energy- 
resolved photoelectron spectra of hydrogen atom using 
few-cycle laser pulses for both homogeneous and inho- 
mogeneous fields, considering tunneling and multiphoton 
regimes. Furthermore, we solve the classical equations 
of motion of an electron in an oscillating inhomogeneous 
electric field to support our quantum mechanical method. 
Finally, in Sec. IV, we conclude our contributions with a 
short summary and outlook. 

II. THEORETICAL APPROACH 

To study the properties of the ATI phenomenon driven 
by spatial nonhomogeneous fields, we solve the three di- 
mensional Time Dependent Schrodinger Equation (3D- 
TDSE) in the length gauge. The electron momentum dis- 
tribution and energy-resolved photoelectron spectra of an 
atom are calculated from the time propagated electronic 
wave function. 

Our calculations are based on spherical harmonics ex- 
pansion, yj m , considering only the m = terms due to 
the cylindrical symmetry of the problem. The Crank- 
Nicholson method, which is based on a splitting of the 
time-evolution operator that preserves the norm of the 
wave function, is used as the numerical technique. We 
consider the field to be linearly polarized along the z 
axis and the variation of the electric field is linear with 
respect to the position. As a result, the coupling V/(r, t) 
between the atomic system and the electromagnetic ra- 
diation reads 

Vi(r,t) = J dr f • E(r', t) = E z(l + f3z)f(t) sm(ut + 0), 

(1) 

where Eq, uo and <p are the laser electric field amplitude, 
the central frequency and characterizes the carrier enve- 
lope phase (CEP), respectively. The parameter f3 defines 
the 'strength' of the inhomogeneity and has units of in- 
verse length (see also [361-38] ) . For modeling the short 
laser pulses in Eq. (pQ), we use a sin-squared envelope f(t) 
of the form 

/<«>—•(£). P) 

where n p is the total number of optical cycles. As a re- 
sult, the total duration of the laser pulse will be T p = n p r 



3 



where r = 2tt/uj is the laser period. We also assume that 
before switch on of the laser (t = — oo) the target atom 
(hydrogen) is in its ground state (Is), whose analytic 
form can be found in standard textbooks. Within the 
single active electron approximation, however, our nu- 
merical scheme is tunable to treat any complex atom by 
choosing the adequate effective (Hartree-Fock) potential, 
and finding the ground state by the means of numerical 
diagonalization. 

The ATI spectrum is calculated using the time depen- 
dent wave function method developed by Schafer and Ku- 
lander (see [48| for more details). As a preliminary test, 
for ensuring the consistence of our numerical simulations, 
we have checked out our calculations with the results pre- 
viously obtained in Ref. [48|. The comparison confirms 
the high degree of accuracy of our calculations as shown 
in Fig. [TJ 
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FIG. 1. (Color online) Photoelectron spectrum resulting from 
our 3D TDSE simulations (in red) and superimposed (in 
black) with the ATI results calculated by Schafer and Ku- 
lander in Ref. [48]. The laser wavelength is A = 532 nm and 
the intensity is / = 2 x 10 13 W/cm 2 (see Fig.l in [48[ for more 
details. The superimposed plot has been extracted from Fig. 
1 of this cited reference. 

For calculating the energy-resolved photoelectron spec- 
tra P(E) and two-dimensional electron distributions 
T-L(p,0) we use the window function approach developed 
by Schafer [48, 49]. This tool has been widely used, both 
to calculate angle-resolved and energy-resolved photo- 
electron spectra |50j, and it represents a step forward 
with respect to the usual projection methods. 

III. RESULTS 

In this section, we calculate both energy-resolved pho- 
toelectron spectra P(E) and two-dimensional electron 
momentum distributions in order to investigate the in- 
fluence of the inhomogeneities of the field and the sensi- 
tivity of these two measurable quantities to the different 
laser parameters, especially to the carrier envelope phase 
(CEP). The investigations are carried out for both the 
tunneling regime, for which the Keldysh parameters is 



7 ^ 1 (7 = V^p/^Up, where U p = //4cj 2 is the pon- 
deromotive energy and I p the ionization potential), and 
multiphoton regime, for which the Keldysh parameters 
is 7 >> 1. Furthermore, we confirm how in the tun- 
neling regime the CEP, joint with the spatial nonhomo- 
geneities, modify in a particular way both the energy- 
resolved photoelectron spectra and the two-dimensional 
electron momentum distributions as we have shown in 
our previous contribution [40]. On the other hand, we 
show that in the multiphoton regime (7 >> 1), the spa- 
tial nonhomogeneous character of the laser electric field 
hardly affects the analyzed quantities. We also want to 
point out, however, that the frontier between the tunnel 
and multiphoton reg imes appears to be a controversial 
and diffuse issue [5l|, [52| . 



A. Tunneling regime 

We commence by investigating the tunneling regime. 
For this case, we employ a four-cycle (total duration 10 
fs) sin-squared laser pulse with wavelength A = 800 nm 
and two different intensities, namely / = 1.140 x 10 14 
W/cm 2 and I = 5.0544 x 10 14 W/cm 2 . These two inten- 
sities give values for the laser electric field of Eq = 0.057 
a.u. and Eq = 0.12 a.u., respectively. For all the cases we 
chose four different values for the parameter that charac- 
terizes the inhomogeneity strength, namely, f3 = (ho- 
mogeneous case), 0.002, 0.003 and 0.005. In addition we 
also vary the carrier envelope phase <p in Eq. (pp), tak- 
ing (j> = 0, <j) = 7r/2, (j) = 7r and <p = 3tt/2. For all 
the above mentioned cases, we calculate the the energy- 
resolved photoelectron spectra. The results are shown in 
Figures 2 and 3 for I = 1.140 x 10 14 W/cm 2 (7 = 1) and 
I = 5.0544 x 10 14 W/cm 2 (7 = 0.475), respectively. 
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FIG. 2. (Color online) Energy-resolved photoelectron spectra 
P(E) calculated using the 3D-TDSE for an hydrogen atom 
(I p — —0.5 a.u.). The laser parameters are / = 1.140 x 10 14 
W/cm 2 (E = 0.057 a.u.) and A = 800 nm. We have used a 
sin-squared shaped pulse with a total duration of four optical 
cycles (10 fs). (a) ft = (homogeneous case), (b) f5 — 0.002, 
(c) p = 0.003 and (d) f3 = 0.005. In all the panels, black line: 
= 0; blue line (j) = 7r/2; green line: (j) = ty and red line: 
= 3tt/2. 

For the homogeneous case, the spectra exhibits the 
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FIG. 3. (Color online) Idem Fig. 2 but for a laser intensity 
I = 5.404 x 10 14 W/cm 2 (E = 0.12 a.u.). 



usual distinct behavior, namely, the 2U P cutoff (« 13.6 
eV in Fig. 2(a) and w 57.4 eV in Fig. 3(a)) and the 
10U P cutoff (« 68 eV in Fig. 2(a) and w 300 eV in 
Fig. 3(a)). The former cutoff corresponds to those elec- 
trons that, once ionized, never return to the atomic core 
(the so-called direct electrons), while the latter one cor- 
responds to the electrons that, once ionized, return to 
the core and elastically rescatter (the so-called rescat- 
tered electrons). Classically, it is a well established ar- 
guments that the maximum kinetic energies Ek of the 
direct and the rescattered electrons are Et 
E r 



- 2U P and 

10U P , respectively (see below for more details). 
In a quantum mechanical approach, however, it is possi- 
ble to find electrons with energies beyond the lOUp cut- 
off, although their yield drops several orders of magnitude 
depending strongly on the atomic species studied [l|. 




FIG. 4. (Color online) Two-dimensional electron momen- 
tum distributions (logarithmic scale) in cylindrical coordi- 
nates (k Z jkr) using the exact 3D-TDSE calculation for an 
hydrogen atom. The laser parameters are I = 1.140 x 10 14 
W/cm 2 (E = 0.057 a.u.) and A = 800 nm. We have used a 
sin-squared shaped pulse with a total duration of four optical 
cycles (10 fs) with = 0. (a) f3 = (homogeneous case), (b) 
P = 0.002, (c) p = 0.003 and (d) = 0.005. 

Experimentally speaking, both the direct and rescat- 
tered electrons contribute to the energy-resolved photo- 



electron spectra. It means for tackling this problem both 
physical mechanisms should be included in theoretical 
model. In that sense, the 3D-TDSE, which can be consid- 
ered as an exact approach to the ATI problem for atoms 
in the single active electron approximation, is the ad- 
equate tool to predict the P(E) in the whole range of 
electron energies. On the other hand, the most energetic 
electrons, i.e., those with kinetic energies Ek ^> 2U P , 
are commonly used to characterize the CEP of few-cycle 
pulses. Consequently, a correct description of the elec- 
tron rescattering mechanism is needed. 




FIG. 5. (Color online) Idem Fig. 4 but for = tt/2 

For the spatially nonhomogeneous cases, the positions 
of the direct and the rescattered electron cutoffs are 
extended towards larger energies. For the rescattered 
electrons, this extension is very noticeable. In fact for 
f3 = 0.005 with Eo = 0.12 a.u., it reaches values close 
to « 700 eV [Fig. 3(d)] against the w 300 eV shown 
by the homogeneous case. Another new feature present 
for all the nonhomogeneous cases is the strong sensitivity 
of the P(E) to the carrier envelope phase (CEP). This 
behavior can be clearly noticed by comparing the panels 
(a) of Figs. 2 and 3 (i.e. the homogeneous case) with 
the rest of the plots. It is clear that for the homogeneous 
case only two curves are present, due to the fact that the 
P(E) for (j) = and (<p = tt/2) are identical to <p = tt 
and (cj) = 37r/2), respectively. On the other hand, for all 
the nonhomogeneous cases it is possible to clearly dis- 
tinguish the 4 cases, i.e. = 0, </> = 7r/2, </> = 7r and 
4> = 37r/2. Indeed, this particular characteristic of the 
P(E) for nonhomogeneous fields could make them a new 
and better CEP characterization tool. 

It should be noted, however, that other well-known 
and established CEP characterization tools, such as, 
for instance, the forward-backward asymmetry or two- 
dimensional electron momentum distributions should 
complement the P(E) measurements [l|. Furthermore, 
the utilization of nonhomogeneous fields would open the 
avenue for the production of high-energy electrons, reach- 
ing the keV regime, if a reliable control of the spatial 
and temporal shape of the laser electric field is attained. 
Investigations in such direction has already started (see 
e.g. [29j and references therein). 

A deep analysis of the electron distributions for atomic 
ionization produced by an external laser field can be per- 
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FIG. 6. (Color online) Idem Fig. 4 but for (j) — tt 



formed in terms of the two-dimensional electron momen- 
tum distribution. The exact solution of the three dimen- 
sional Schrodinger equation (3D-TDSE) provide us with 
an excellent tool to analyze in detail how the two compet- 
ing fields, namely the laser electric field and the Coulomb 
potential, modify the electron wavepacket of the released 
electron. In Figs. 4-7 we calculate two-dimensional elec- 
tron momentum distribution for a laser field with an in- 
tensity of I = 1.140 x 10 14 W/cm 2 (E = 0.057 a.u), 
A = 800 nm and different values of the the (3 parameter: 
panel (a) (3 = (homogeneous case); panel (b) (3 = 0.002; 
panel (c) (3 = 0.003 and panel (d) (3 = 0.005. We employ 
a few-cycle laser pulse with 4 total cycles (10 fs) and vari- 
ous values of the carrier envelope phase (CEP) parameter 
(j). Figs. 4-7 show the cases with <\> = 0, <p = 7r/2, <p = tt 
and (j) = 37r/2, respectively. 
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FIG. 7. (Color online) Idem Fig. 4 but for = 3tt/2 

Here, we concentrate our analysis on the low energy 
region of the distributions in order to study how the in- 
homogeneities of the laser electric field affect the angular 
electron yield. This region shows the usual bouquet-type 
structures (see [ll[ for details) with noticeable modifica- 
tions for the nonhomogeneous cases. 

Furthermore, the low-energy electrons appear to be 
strongly influenced by the spatial inhomogeneity of the 
laser electric field (see panels (b)-(d) of Figs. 4-7). We 
also can observe how the bouquet structures present in 
the homogeneous case disappear for particular values of 
4> (see e.g. Fig. 6(d)). 

In order to complete our investigations, we calculate 
two-dimensional electron momentum distributions by in- 
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FIG. 8. (Color online) Two-dimensional electron momen- 
tum distributions (logarithmic scale) in cylindrical coordi- 
nates (k z ,k r ) using the exact 3D-TDSE calculation for an 
hydrogen atom. The laser parameters are / = 5.0544 x 10 14 
W/cm 2 (E = 0.12 a.u.) and A = 800 nm. We have used a 
sin-squared shaped pulse with a total duration of four optical 
cycles (10 fs) with = 0. (a) (3 = (homogeneous case), (b) 
(3 = 0.002, (c) (3 = 0.003 and (d) = 0.005. 



creasing the laser field intensity to / = 5.0544 x 10 14 
W/cm 2 (Eo = 0.12 a.u). The results are depicted in 
Figs. 8-10 for <j> = 0, (j) = tt/2, (j) = tt and <j> = 3tt/2, 
respectively. Here, panel (a), (b), (c) and (d) represent 
the cases with (3 = (homogeneous case), (3 = 0.002, 
(3 = 0.003 and (3 = 0.005, respectively. 




FIG. 9. (Color online) Idem Fig. 8 but for = tt/2 




FIG. 10. (Color online) Idem Fig. 8 but for = tt 



Following the trend observed in the previous studied 
case, we see strong modifications produced by the spa- 
tial inhomogeneities in both the angular and low-energy 
structures. Here, the modification seems to be even more 
pronounced. For instance, the yield for electrons with 
k z < -1 a.u. for cj) = (Fig. 8(a)) and (j) = tt/2 
(Fig. 9(a)) drops several orders of magnitude. The sig- 
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nificant drop of the yield of the electron emission between 
k z < — 1 a.u and k z > a.u. opens a new approach to 
characterize the CEP. 

We now employ a classical model in order to explain 
and characterize the extension of the energy-resolved 
photoelectron spectra. According to the simple-man's 
model [53j| the physical mechanism behind ATI process 
will be as follows: at a given time, that we call ionization 
time ti, an atomic electron is released or born at the po- 
sition z = with zero velocity, i.e., z(tj) = 0. This elec- 
tron now moves only under the influence of the oscillat- 
ing laser electric field (this model neglects the Coulomb 
interaction with the remaining ion) and will reach the 
detector either directly or through the process known as 
rescattering. By using the classical equation of motion 
of an electron moving in an oscillating electric field, it is 
possible to calculate the maximum energy of the electron 
for both the direct and rescattered processes. 
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FIG. 11. (Color online) Idem Fig. 8 but for = 3tt/2 

The Newton equation of motion for the atomic elec- 
tron in the laser electric field can be written, using the 
functional form of Eq. (pQ) , as follows 



z(t) = -VzVi(r,t) 

= E(t)[l + 2Pz(t)] 



(3) 
(4) 



where we have collected the time-dependent part of the 
electric field in E(i), i.e., E(t) = E f(t) sm(ut + (j)). 

In the limit where (3 = in Eq. (|3j), we recover the 
homogeneous case. For the direct laser ionization, the 
kinetic energy of an electron released or born at time ti 



E d k = 



[z( ti )-z(t f )Y 



(5) 



where tf is the end time of the laser pulse. For the rescat- 
tered laser-ionized electron, in which the electron returns 
to the core at a time t r and reverses its direction, the ki- 
netic energy of the electron yields 



El 



[z(U) + z(t f )-2z(t r )Y 



(6) 



For homogeneous fields, Equations (j5]) and (j6j) be- 
come the usual expressions E% 

[A(tj)+A(t f )-2A(t r )] 2 



-, with A(t) being the laser vector po- 
tential A(t) = — J 1 E(t')dt' ', respectively. For the case 
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FIG. 12. (Color online) Numerical solutions of the Newton 
equation [Eq. ([3])] plotted in terms of the direct and rescat- 



and El 



respectively. 

2 



The 



tered electron kinetic energy, E k 

laser parameters are I = 1.140 x 10 14 W/cm 2 (E = 0.057 
a.u), A = 800 nm and = 0. We employ a few-cycle laser 
pulse with 4 total cycles (10 fs). Different panels correspond 
to various values of (3 (see labels). Green filled circles: direct 
electrons; red filled circles: rescattered electrons. 



with /3 = 0, it can be shown that the maximum value for 
E% is 2U P while for E r k it is 10U P Q. These two values 
appear as cutoffs in the energy resolved photoelectron 
spectrum as can be observed in panels (a) of Figs. 2 and 
3. 

In Figs. 12-19, we present the numerical solutions of 
Eq. (|3]), which is plotted in terms of the kinetic energy 
(in eV) of the direct and rescattered electrons. 

Figures 12-15 are for a laser intensity of / = 1.140 x 
10 14 W/cm 2 (E = 0.057 a.u.) meanwhile in Figs. 16-19 
the laser intensity is I = 5.404 x 10 14 W/cm 2 (E = 0.12 
a.u.). Figures 12 (16), 13 (17), 14 (18) and 15 (19) are for 
(j) = 0, (j) = 7r/2, (j) = 7r, (j) = 37r/2, respectively and for 
different values of the f3 parameter ((3 = (homogeneous 
case), (3 = 0.002, (3 = 0.003 and (3 = 0.005, from top to 
bottom). 

From the curves for (3 7^ we can observe the strong 
modifications that the nonhomogeneous character of the 
laser electric field produces in the kinetic energy of the 
laser-ionized electron. Furthermore, it is possible to ob- 
serve how the kinetic energy is more sensitive to the CEP. 
If we analyze, for instance, the top panels of Figs. 12-15 
(i.e the homogeneous case) we conclude that the shape 
of both the kinetic energy of the direct and rescattered 
electron are identical for (j) = (0 = tt/2) and </> = tt 
(0 = 37r/2). On the other hand, for the different values 
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FIG. 13. (Color online) Idem Fig. 12 for = tt/2. 
E =0.057 a.u. 4>=tu 
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FIG. 14. (Color online) Idem Fig. 12 for = tt. 



of f3 the kinetic energy has a unique shape for a given 
value of <p. 

The particular features present for f3 ^ are related to 
the changes in the laser-ionized electron trajectories (for 
details see e.g. [37-39]). In summary, the electron trajec- 
tories are modified in such a way that now the electron 
ionizes at an earlier time and recombines later, and in 
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FIG. 15. (Color online) Idem Fig. 12 for = 3tt/2. 



this way it spends more time in the continuum acquir- 
ing energy from the laser electric field. Consequently, 
higher values of the kinetic energy are attained. This 
distinct behavior is more evident for Eo = 0.12 a.u. and 
/3 = 0.005, but it appears to some degree for all the stud- 
ied cases. 

A similar behavior was observed recently in above 
threshold photoemission (ATP) using metal nanotips. 
According to the model developed in Ref. [3l| the strong 
localized fields modify the electron motion in such a way 
to allow sub-cycle dynamics. 

In our approach, however, we include the full picture 
of the ATI phenomenon, namely both direct and rescat- 
tered electrons are considered (in Ref. [3l| only direct 
electrons are taken into account) and consequently the 
characterization of the dynamics of the photoelectrons is 
more complex. Nevertheless, the higher kinetic energy 
of the rescattered electrons is a clear consequence of the 
strong modifications the laser electric field produces in 
the region where the electron dynamics takes place, as in 
the above mentioned case of ATP. 



B. Multiphoton regime 

For the multiphoton case, we consider a few-cycle laser 
pulse with 6 complete optical cycles and Eq = 0.05 a.u. 
(I = 8.775 x 10 13 W/cm 2 ) and u = 0.25 a.u. (A = 182.5 
nm). In here, the Keldysh parameter is 7 = 5, indicating 
the predominance of the multiphoton process [11]. We 
have computed the P(E), two-dimensional electron dis- 
tributions and the classical electron energies for all the set 
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E Q =0.12a.u. <H) 




Time (optical cycles) 

FIG. 16. (Color online) Numerical solutions of the Newton 
equation [Eq. [3] plotted in terms of the direct and rescattered 
electron kinetic energy, E% and El, respectively. The laser 
parameters are I = 5.0544 x 10 14 W/cm 2 (E = 0.12 a.u), 
A = 800 nm and = 0. We employ a few- cycle laser pulse 
with 4 total cycles (10 fs). Different panels correspond to 
various values of f5 (see labels). Green filled circles: direct 
electrons; red filled circles: rescattered electrons. 



of cases presented in Sec. III. A. In this paper, however, 
we just present the most extreme case with an inhomo- 
geneity factor of f3 = 0.005 and CEP of qb = tt/2. These 
results are presented in Figs. 20, 21 and 22 for P(E), 
two-dimensional electron distributions and the classical 
electron energies, respectively. The P(E) exhibits the 
usual multiphoton peaks 0, [48| and the inhomogeneity 
of the field does not play any significant role. In the 
whole range, the values of the yields have a difference of 
less than 5% and in a logarithmic scale this is hard to 
discern. 

The two-dimensional electron distributions are also the 
same in terms of shape and magnitude for both homoge- 
neous and inhomogeneous cases, as shown in Fig. 21. It 
means the differences introduced by the spatial nonho- 
mogeneous character are practically imperceptible. We 
should note that our calculation is practically identical 
to the one presented in [ll|. 

The numerical solutions of Eq. (J3j) as function of the ki- 
netic energy (in eV) of the direct and rescattered electron 
is depicted in Fig. 22. In here, we could also observe, in 
support to our quantum mechanical calculations, that the 
inhomogeneity of the field does not change the electron 
energies in both the direct and rescattered processes. 

In general, we do not find noticeable differences in 
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FIG. 17. (Color online) Idem Fig. 16 for = tt/2. 
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FIG. 18. (Color online) Idem Fig. 16 for = tt. 



these observable quantities for both variations in the 
CEP and the strength of the inhomogeneity parameter 
f3. As a result, we conclude that in the multiphoton 
regime the modifications introduced by the spatial in- 
homogeneities field do not produce appreciable modifica- 
tions in the electron dynamics and consequently in the 
measurable quantities. In addition, the laser-ionized elec- 
tron, in both the direct and rescattered processes, has a 
very small kinetic energy in the multiphoton regime due 
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E=0.12 a.u. (J)=3ti/2 
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FIG. 19. (Color online) Idem Fig. 16 for = 3tt/2. 
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FIG. 20. (Color online) Energy-resolved photoelectron spec- 
tra P(E) calculated using the 3D-TDSE for an hydrogen atom 
(I p = —0.5 a.u.). The laser parameters are Eo — 0.05 a.u. 
(I = 8.775 x 10 13 W/cm 2 ) and oj = 0.25 a.u. (A = 182.5 nm). 
We have used a sin-squared shaped pulse with a total dura- 
tion of six optical cycles and <j> — tt/2. Red line: homogeneous 
case (J3 = 0); blue line: f3 = 0.005;. 




k z (a.u.) 



FIG. 21. (Color online) Two-dimensional electron momen- 
tum distributions (logarithmic scale) in cylindrical coordi- 
nates (k z ,k r ) using the exact 3D-TDSE calculation for an 
hydrogen atom. The laser parameters are Eo — 0.05 a.u. 
(J = 8.775 x 10 13 W/cm 2 ), u = 0.25 a.u. (A = 182.5 nm) and 
cj) — 7r/2. We employ a laser pulse with 6 total cycles. Panel 

(a) corresponds to the homogeneous case (f3 = 0) and panel 

(b) is for P = 0.005. 



IV. CONCLUSIONS AND OUTLOOK 

We have extended our previous studies of ATI pro- 
duced by nonhomogeneous fields using the three dimen- 
sional solutions of the TDSE. We have modified the 3D- 
TDSE to model the ATI phenomenon driven by spatial 
nonhomogeneous fields by including an additional term 
in the laser-atom coupling. In the tunneling regime we 
predict an extension in the cutoff position and an in- 
crease of the yield of the energy-resolved photoelectron 
spectra in certain regions. In addition, both the pho- 
toelectron spectra and the two-dimensional electron mo- 
mentum distributions appear to be more sensitive to the 
carrier envelope phase of the laser electric field. This 
feature indicates that the ATI produced by spatial in- 
homogeneous field could be a good candidate for few- 
cycle laser pulse characterization. Furthermore, our pre- 
dictions pave the way for the production of high-energy 
photoelectrons, reaching the keV regime, using plasmon 
enhanced fields. In the multiphoton regime, we show that 
both the P(E) and the two-dimensional electron distribu- 
tions are hardly affected by the spatial nonhomogeneities 
of the laser electric field. Our quantum mechanical cal- 
culations are supported by the classical simulations. In 
particular, the P(E) characteristics are reasonably well 
reproduced by simulations based on classical physics. 



to low intensity field. Indeed, in this case the maximum 
energy after the rescattering process has taken place is 
w 3 eV. As a result, it is reasonable to have a very small 
or almost no differences between the final kinetic ener- 
gies, when a spatial inhomogeneity of small strength is 
present. 
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FIG. 22. (Color online) Numerical solutions of the Newton 
equation [Eq. (|3])] plotted in terms of the direct and rescat- 
tered electron kinetic energy, E% and El, respectively. The 
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employ a laser pulse with 6 total cycles. Panel (a) corre- 
sponds to the homogeneous case (f3 = 0) and panel (b) is for 
f3 = 0.005. Green filled circles: direct electrons; red filled 
circles: rescattered electrons. 
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